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Abstract. The question of robustness of a network under random "attacks" is treated 
in the framework of critical phenomena. The persistence of spontaneous magnetization 
of a ferromagnetic system to the random inclusion of antiferromagnetic interactions is 
investigated. After examing the static properties of the quenched version (in respect 
to the random antiferromagnetic interactions) of the model, the persistence of the 
magnetization is analysed also in the annealed approximation, and the difference in 
the results are discussed. 
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1. Introduction 

The investigation on resilience of networks, which examines the persistence/breakdown 
of some global properties of a graph under, for instance, removal of vertices or edges, is 
known to have practical importance. Many "real networks" (Internet, highways, many 
biological systems, et cetera) depend on the fact that there exist links between the nodes 
to ensure their functionality. A damage that breaks the interconnection between vertices 
can trigger a profound impact on the network]!]. The connection of these questions to 
the percolation theory is clear, which has been a powerful tool to tackle these problems [2] 
(see also 011]). 

In this work, the question of resilience of networks is analysed in the framework 
of critical phenomena of magnetic systems, and is not directly related to the classical 
problems concerning the existence of giant components in a graphp2,[6]. 

Critical phenomena have been exhaustively exploited in networks]?! [8], and it 
was immediately noticed that it displays a quite rich thermodynamical behaviour 
when compared to regular lattice structures [TTJ [121 H31 EH] or gives place to nontrivial 
effects [T5l [16]. This shows the role of the topology of the underlying graph on critical 
behaviour. Moreover, networks constitute a path toward a more "realistic" system of 
mean - field models, the main cause being the introduction of (finite) connectivity as a 
parameter pH M, CEHJ M> EH E2] . 

This work will focus on the persistence of a order imposed by a background 
enviroment to random inclusion of opposite effects on the network. More precisely, 
starting from a ferromagnetic system, some of its interactions are weakened (or even 
changed) by antiferromagnetic ones. The criterion for the "functionality" of the network 
will be the ferromagnetic order of the model, and by "resilience" it means how strong 
is the spontaneous magnetization against random introduction of antiferromagnetic 
interaction that contributes to the disorder. The model is detailed in section II, and 
the thermodynamical analysis of its quenched version is given in section III, followed 
by a characterization of the order - disorder critical line in the annealed approximation 
in section IV. The critical line that determines the breakdown of the ordered phase is 
discussed in section V, and the last section is devoted to conclusions. 



Consider a (finite) graph Ajv (| Ajv| — N) where each vertex i 6 Ajy allocates an Ising 
spin a x G { — 1,1}. The Hamiltonian is given by 



where the first term represents the ferromagnetic interaction (J/N > 0) between every 
pair of spins and the second sum is the antiferromagnetic interaction (—J a < 0) between 
some pair of spins. Therefore, the network is a complete graph where the edges represent 



2. Model 




(1) 



x<y 
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interactions. Note that the antiferromagnetic interactions do not constitute new edges 
on the graph themselves, but they weaken the pre - existent ferromagnetic edges 
competitively. The matrix {a xy } decides if two vertices, x and y, have antiferromagnetic 
interaction (a xy = 1) or not (a xy = 0), and the couplings {a xy } are independent and 
identically distributed random variables that obey the distribution 

V{a xy ) =(l-^j 5 axyfi + 1^,1 • (2) 

It is known that the quenched version of model (PQ) without the ferromagnetic 
term (the first sum) is a frustrated system with a spin glass phase at low 
temperatures [23 ] (since the model has a single sublattice). Therefore, between the two 
effects of antiferromagnetic interactions, namely the contribution to disorder and a 
contribution to ferromagnetic order (due to the effect of frustration arisen in some loops 
of the graph), the former dominates over the later. 

As usual, the partition function is 

Z(P, {a xy }):= Tr K} e-«^»-» , (3) 

where (3 is the inverse of the temperature T (ks = 1) and the trace indicates the 
sum over 2 N states. This partition function also depends on the configuration of the 
antiferromagnetic interactions. Finally, given a function g = g{{a xy }) of the set of 
random variables {a xy }, the average of g over the configuration {a xy } will be denoted 

as 

(g) ■= / g({a xy })Y[V({a xy })da xy . (4) 

x<y 

3. Thermodynamics 

The (quenched) free energy f q of the model is evaluated through the replica trick 

fM = - lim ^QnZ(Jf)) = -lim lim ^ln(i? n (/?)>, (5) 

N^oo ply n-*0N^oo pl\n 

where besides the analytic continuation in n G N — > n G R, the order of the limits 
iV — >• oo and n — > was changed as usual. Then, one should compute {Z n (/?)), which 
can be casted as 



(f)"nH^[-f^ 

\a-l R ) a-l 

n 



(6) 



Throughout this work, upper and lower indices at spin variables indicate replica index 
(greek letter) and site position (roman letter), respectively. 
Introducing the order parameter [ 



^ := ^E h^},{^} > ( 7 ) 



N 

xGAjv 



Diluted antiferromagnet in a ferromagnetic enviroment 



4 



where ^{a^},{^ a } '■= Yla=i the equation §6§ can be written as 



(^) f (n-)/^ 



-N<t> q [i>,iP](i3,{\°'}) 



(8) 



where 



$)(p, {A*}) := P - + ^-it (A°) 2 + Tr^Mfo) - 



a=l 



P 



and 



C[^](/3,{A Q }):=Tr {(TQ} exp 



a=l 



(9) 



(10) 



Apart from a factor /3, ^["0, V'K/^, {A a }) is just the variational free energy. The 
equation (jSJ) suggests that one should invoke the saddle - point method to determine 
the stationary free energy. The extremum conditions, necessary to ensure the infimum 
of (f) q over the suitable functions (ip and 0) and variables ({A a }), are 



0(/i) 



1 



exp 



a=l 



ip(fi) = p Tr {r a } ^(r) exp 



a=l 



1 



Ti{ a aya a exp 



a=l 



CM CM**}) 

In an attempt to solve the above equations, one should cast the replica symmetric 
Ansatz 



[2 cosh(/3/i)f 



W = i>(I2 n a=i» a )=P dyQ(y)-, 



(12) 



[2cosh(/fy)] r 



A a = A, Vae{l,---,n} 
where P and Q are probability distributions. 
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A 



oo _ / r 
P 



r=0 ' \j=l 



e p Hv\ I ny i <>x i > 

,=1 R 

cosh (/3Ja — 



/3x ( JA + 



fc=i 



(13) 



cosh ((3J A + /3h)_ 

^Etf ny^w tanh 



r=0 \ j=l 

the (quenched) free energy f g is evaluated as 

MUD PJy ' P 



k=l 



+ |ln2 + p J dhP(h) J dyQ(y) In cosh \J3 (h + y)) - 

R R 

-| J dhxPihx) J dh 2 P{h 2 ) In je^ 1 cosh [/3 (/i 2 - J A )\ + e"^ 1 cosh [/3 (h 2 + J A 



-ln2-e" p ^ 



r=0 



V 



s=l 



1 [ / dy s Q(y s 



In cosh 



/?( JA + J> 



(14) 



the extremum condition leads P and Q to satisfy 



P(h) = e-pJ2^[Ul d Vj Q{y j )\8[h 



r=0 ' \j=l 



k=l 



(15) 



Q( y ) = / dhP(h)5 (v+jj tanh" 1 [tanh(/3J A ) tanh(/3/i)] J 



and A is calculated through the last equation of f|T3|) . The above set of equations ffT5l) 
can be unified as 



r=0 ' \j=l 



xdU- 



X 



1 r 

JA — — tanh" 1 [tanh(/3Jyi) tanh(/3/ifc)] 



k=l 



(16) 



In the context of replica method, the formulas for the magnetization, m, and the 
spin - glass order parameter, q, are given by 



m 



1 1 " / n 

lim — Y lim - V Trr^n^ ( TT e" 



-/^({o-ZK-K*}) 



TV— >oo TV ^— ' n— >o n 



(17) 



\7=1 



and 



xSAjv V 7 a^6> \7=1 / 
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Using (TT5|) . (TTTT) and (IT5|) . it is possible to see that A = m, 

m= I dhP(h)timh(ph) and q = / dhP{h) tanh 2 (/3/i) , (19) 



as usual. The critical line is determined in the neighborhood of m ~ and g ~ 0. In 
this regime, the field /i is expected to be narrowly distributed around h = 0. Then, the 

Ansatz 

i* := / dhP(h)h k = 0{e k ) , |e|<l, (20) 



is introduced in the equation ffl6l) to evaluate the transition lines [13J. This Ansatz leads 
to m = 0(e) and q = 0(e 2 ), which means that the line of transition from ferromagnetic 
to disordered phase is governed by 0(e) in the equation f[T5|) and the transition from 
paramagnetic to spin - glass phase by the order C(e 2 ) (e = is assumed in this case). 
As a result of these calculations, one has 

j3 c J = 1 + ptanh(/? c Jyi) (order - disorder transition line) , (21) 

and the transition line from spin - glass to paramagnetic phase is evaluated as 

=7777^1* P>1- (22) 



The phase diagram of the model, generated by the equations ( 12TT) and ( 1221) . is 
presented in FIG 1. The "m 7^ 0" phase displays spontaneous magnetization, and it 
may be a combination of a ferromagnetic and a mixed phases. The exact scenario of 
the "m 7^ 0" phase can be established from a stability analysis (25^, which will not be 
provided here, since hereafter this work will focus on the critical line that separates the 
ordered phase and disordered one in comparison with the annealed version of the model. 

4. Annealed average 

In the annealed approximation, the free energy f a is written as 

/.(# = - Hm -±-ln(Z((3)). (23) 

N-*oo ply 

The problem can be solved in the annealed approximation in a simpler way. Although 
the interesting case is the quenched one (which was provided in the previous section), 
the results will be derived for comparison. 
It is straighforward to show that 

(Z(/?)> = Tr W} exp I - ^ + ^ cosh(/3J A ) + ^ [(3J -psmh((3J A )] ( £ + 
+(9(1)1, (24) 



X To be more precise, the stability analysis can also change the exact location of the line that separates 
the spin - glass phase and the paramagnetic phase; however, the phase diagram showed in FIG 1 is 
believed to be qualitatively correct. 
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Figure 1. Phase diagram with J a/ J = 1- 



and introducing the order parameter 



m N :-- 



N ? ° x 



it is easy to show that 

(Z((3)) ~ / dm N e- N(t>a ^ mN) 



for sufficiently large N, where 



<t> a (/3,m N ) := - 



p pcosh(/3J^) /3J 2 psinh(/3J,4 



(25) 



(26) 



m N + -^-ln 



■itln ( 1 + rriN 



1 — 772 TV 



+ 



+-ln(l-m^) -ln2. 



(27) 



Therefore, in the thermodynamic limit, with m := limjv^oo mjv, the free energy can 
be written as 

/ a (/3) = ^inf{0 o (/3,m)}, (28) 

p m 

where the infimum of <fr a is achieved from the solutions of the extremum condition 

m = tanh (3Jm — psinh(/9 Jj\)m . (29) 

This equation allows one to obtain the critical temperature j3 c , which then obeys 

(3 C J= l+psinh(&J A ), (30) 

and is cleary different from the correspondent expression (l2Tj) from the quenched 
situation. 
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5. Breakdown of the spontaneous magnetization 

This section will analyse the different behaviour of the critical line evaluated in the 
quenched and annealed approaches, which are 

P c J = 1 + ptanh(/? c J J 4) (quenched case) (31) 

and 

(3 C J = 1 + psmh(j3 c JA) (annealed case) . (32) 

Firstly, it is easy to see that if p = 0, one recovers the ferromagnetic mean - field 
result (3 C J = 1 in both cases, as it should be. Now, let p assume nonzero values. Actually, 
given p and Ja/J, the equation (132]) may yield two roots, but only the physically 
reasonable one for the order - disorder transition is chosen. The numerical solutions 
of the equations are plotted in FIG 2 as a function of p. As one can see in the figure, 
for sufficiently small J a (in the sense that (3 c Ja "C 1), the critical temperature describes 
a line T c j J ~ 1 — p(Ja/ J) in both cases. 
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Figure 2. Dependence of the critical temperature on the mean connectivity p (left: 
quenched case; right: annealed case). 

However, for a fixed value of J a/ J, one can see the difference from the quenched and 
annealed cases as the mean connectivity p increases. The spontaneous magnetization is 
broken for sufficiently large values of p in an abrupt way in the annealed case. On 
the other hand, the critical temperature for the quenched case decays slowly, and 
reaches T c = for p — > oo only (see ( f3~TI) ). One should remember, initially, that each 
vertex links to another N — 1 with ferromagnetic interaction J/N = 0(N~ l ) and an 
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antiferromagnetic interaction (which is much stronger in the sense that J a = 0(1)) with 
probability p/N = 0(N- L ). 

In the quenched case, where the configuration {a xy } remains frozen during an 
observational time, only an infinitesimal fraction (of 0(N~ r )) of the links have 
antiferromagnetic interactions. Suppose that the edge xy (shared between the vertices 
x and y) is one of them; this edge has a ferromagnetic component of intensity J/N and 
an opposite effect of J a (3> J/N). Although the antiferromagnetic part exceeds the 
ferromagnetic one, the solely effect is make just the spins a x and a y having opposite 
signs. Therefore, despite the fact that all the antiferromagnetic interactions dominates 
over the ferromagnetic interactions in the links where a xy = 1, the total number of such 
edges is much smaller than the total number of edges of the network (the thermodynamic 
limit is taken for a fixed value of p), which then becomes predominantely ferromagnetic. 
This is the cause of persistence of spontaneous magnetization for any finite p, as shown 
in FIG 2 (left). 

On the other hand, in the annealing approximation, where the antiferromagnetic 
links fluctuate during the observation time, one sees an "averaged" antiferromagnetic 
interaction between spins (vertices). This means that, although the mean connectivity 
p is fixed, the effective number of edges with an antiferromagnetic interaction is much 
larger (and the effective intensity is also smaller than J a). Heuristically speaking, the 
intensity J a is better distributed over the edges (differently from the quenched case), 
and the antiferromagnetic effect is better exploited in the annealed case, which makes 
the magnetization vanishes even for finite values of p. 

6. Conclusions 

Throughout this work, the static properties of a diluted antiferromagnet on a 
ferromagnetic background was examined, with particular emphasis on the breakdown 
of the spontaneous magnetization of the system. The phase diagram displayed a 
nonzero magnetization at low temperature regime for any finite mean connectivity p of 
antiferromagnetic interactions. The disordered phase is constituted by a paramagnetic 
and a spin - glass phase. 

The critical temperature of the order - disorder transition, which indicates the 
breakdown of the magnetic order of the model, was determined from both the quenched 
and annealed approaches. The main difference, noticed from FIG 2, relies on the 
fact that the spontaneous magnetization vanishes for finite values of p in the annealed 
approximation. This phenomena is observed due to the rapid fluctuation of the random 
variables {a X y}, which distributes more efficiently the antiferromagnetic interactions over 
the whole graph. This means that in the present work, networks are resilient to non - 
fluctuating random "attacks" (even strong ones), while they are weaker to "annealing 
attacks", which turns the antiferromagnetic interaction more accessible to the edges, 
although weakening its mean strength. 
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